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Abstract. We examine the constraints on models of weakly interacting massive particle 
(WIMP) dark matter from the recent observations of the Galactic Center by the High En- 
ergy Spectroscopic System (HESS) telescope. We analyze canonical WIMP annihilation into 
Standard Model particle final states, including bb, ti and W^W~ . The constraints on an- 
nihilation into bb is within an order of magnitude of the thermal cross section at ~3 TeV, 
while the r+r^ channel is within a factor of ~2 of thermal. We also study constraints 
on Sommerfeld-enhanced dark matter annihilation models, and find that the gamma-ray 
observational constraints here rule out all of the parameter space consistent with dark mat- 
ter annihilation interpretations of PAMELA and the Fermi-LAT e~^e~ spectrum, in specific 
classes of models, and strongly constrains these interpretations in other classes. The gamma- 
ray constraints we find are more constraining on these models, in many cases, than current 
relic density, cosmic microwave background, halo shape and naturalness constraints. 
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1 Introduction 

The existence of cosmological dark matter has been well-established by observations of galaxy 
clusters, galaxy rotation curves, the cosmic microwave background (CMB), and large-scale 
cosmological structure. However the identity of the dark matter has remained a fundamental 
unsolved problem in cosmology and particle physics for nearly 80 years [1]. Several particle 
candidates have been proposed that could account for the dark matter (for a review, see, 
e.g. [2]). One well-motivated dark matter candidate is a weakly-interacting massive particle 
(WIMP), which can naturally produce a relic abundance at the observed dark matter density. 
Thermal production of dark matter prefers a scale of the dark matter cross-section at {crp^v) ~ 
3 X 10~^^ cm^ s~^. This annihilation rate into Standard Model particles results in energetic 
gamma-ray production through the hadronization of quarks, bremsstrahlung of leptons, or 
directly into two gammas through higher-order processes. This leads to the so-called method 
of "indirect detection" of dark matter, constraining the dark matter mass, annihilation cross- 
section, and annihilation spectrum through a search for the Standard Model byproducts of 
WIMP annihilation. 

The High Energy Stereoscopic System (HESS) telescope has strong sensitivity to high- 
energy gamma-rays such as those from high-mass WIMPs [3]. HESS consists of an array of 
atmospheric Cerenkov telescopes in Namibia designed to search for high-energy gamma-rays, 
with 960 pixels per telescope at a resolution of 0.16° per pixel [4]. Specifically, HESS is 
sensitive to gamma-ray energies from a few hundred GeV to a few tens of TeV. Previously, 
studies of the Galactic Ridge for \i\ < 0.8° and \b\ < 0.3° by HESS [5] have been used to 
limit the dark matter cross-section at high masses (0.5-30 TeV) [6-13]. However, an even 
more stringent constraint on the dark matter cross-section has been shown to come from a 
new analysis of a region around the Galactic Center (GC) [3]. Here, we analyze in detail the 
constraints arising from this observation. 

The HESS GC analysis uses a refiected background technique to provide a robust back- 
ground region for the gamma-ray signal from the GC. Because the background region is 
further from the GC than the source region, it is expected that the dark matter signal should 
be larger in the source region than the background region, for a dark matter halo profile 
whose density peaks toward the GC. The HESS GC analysis shows no excess gamma-ray 
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signal in the source region over the background region. Therefore, any dark matter annihi- 
lation signal must produce few enough gamma-rays that the gamma-ray flux in the source 
region is indistinguishable from the gamma-ray flux in the background, at the current HESS 
sensitivity. 

There have recently been stacked studies of the gamma-ray signals coming from dwarf 
spheroidal galaxies by the Fermi-LAT which have been used to constrain WIMP annihilation 
cross-sections and, in particular, have excluded thermal low-mass WIMPs in the bb and t~^t~ 
annihilation channels below ~ 30 GeV [14, 15]. Here, we show that the HESS GC analysis 
provides much tighter constraints than the stacked Fermi-LAT dwarf spheroidals for dark 
matter masses above the Fermi-LAT energy window. 

There remains interest in the possibility of dark matter annihilation as the source of the 
excess cosmic ray positron fraction at ~ 10-100 GeV observed by the PAMELA satellite, with 
e~^e~ pairs produced either directly or indirectly in a dark matter particle pair annihilation 
cascade [16-19]. Additionally, features in the higher-energy 10^ to 10^ GeV e+e~ spectrum 
seen by the Fermi-LAT [20] are also consistent with the dark matter annihilation interpre- 
tations of the lower energy positron excess data [10, 21, 22]. A recent study of the e~^e~ 
data from the Fermi-LAT is consistent with the positron excess of the PAMELA satellite 
and shows the spectrum continuing to rise up to at least ~ 200 GeV [23]. 

In order to achieve the dark matter annihilation rate required for these e^e~ signals 
while remaining consistent with the expected thermal production cross section, and to avoid 
an excess in anti-proton observations (which is not seen), the annihilation rate can be en- 
hanced through a low-energy Sommerfeld-enhancement and limited to leptonic modes with 
a light (< 1 GeV) dark-force carrying particle [24-27]. Such an enhanced cross-section from 
a new force is in tension with detailed calculations of the relic abundance of the dark matter, 
so that such a candidate in many cases may not contribute to all of the dark matter [28- 
32]. Such candidates are also constrained by non-thermal distortions of the CMB [29, 33-35] 
and asphericity observed in dark matter halos [30, 31]. Cases of these models remain viable 
given all such constraints [36, 37]. There are also constraints on these models from the ob- 
served diffuse gamma-ray and X-ray backgrounds [38, 39], observations toward the GC by 
Fermi-LAT [68], as well as big-bang nucleosynthesis [40, 41]. 

One model of dark matter with a light dark-force carrying particle and Sommerfeld- 
enhanced cross-section is "eXciting dark matter" (XDM) [24, 42]. XDM was initially pro- 
posed to explain the 511 keV gamma-ray signal from the GC (see ref. [43] for a discussion of 
these signals). In XDM, there are two lowest-energy dark states x X* which have masses 
differing by only a few MeV, with a light gauge bosons mediating excitations from x to 
X* ■ The exchange of many gauge bosons leads to a Sommerfeld-enhanced cross-section much 
larger than that of a thermal relic. The annihilation xx ~^ 4'4' followed by the decay of (p 
into leptons leads to an excess of high-energy electrons and positrons in the GC. This sce- 
nario has also been considered to explain the excess in local positrons seen by the PAMELA 
and Fermi-LAT satellites. Most recently, ref. [36] has interpreted the XDM explanation for 
the PAMELA and Fermi-LAT excesses including constraints from the thermal relic density, 
CMB, self-interaction bounds, and naturalness bounds. Below, we compare the gamma-ray 
constraints from the HESS GC analysis to these other XDM limits. 

The sensitivity of the HESS GC observation to dark matter annihilation, as with all 
Galactic Center observations, depends on the nature of the dark matter density profile. 
Specifically, for the HESS GC background subtraction method, there must be a higher dark 
matter density within the inner <150 pc, in the signal region, than the background subtrac- 
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tion region between approximately 150 pc and 450 pc. As pointed out also in the HESS 
Collaboration analysis [3], if there exists a constant-density core within the inner ~450 pc 
of the Milky Way, no limits on dark matter annihilation can be derived from the HESS GC 
observation since the background subtraction would also remove any equivalent signal. 

There is significant debate in the literature as to the nature of the inner dark matter 
profile of a galaxy such as the Milky Way. Numerical simulations are employed in attempts 
to accurately determine the inner dark matter density profile. The canonically-adopted dark 
matter halo density profile for the case of cold dark matter is the inner-radius scal- 
ing Navarro- Frenk- White (NEW) profile [44]. The highest spatial resolution simulations of 
Milky- Way-type halo formation are pure dark matter halo simulations: Via Lactea H [45], 
GHALO [46] and AQUARIUS (A-1) [47] (which have gravitational softening lengths of 40, 
61 and 20 pc, respectively), which find a peaked density profile down to ~100 pc, with a 
logarithmically-changing slope that is sometimes dubbed an "Einasto" profile. The vast ma- 
jority of studies that include baryons in addition to dark matter have found that baryonic 
effects concentrate and steepen the central dark matter distribution due to adiabatic contrac- 
tion [48, 49], including recent high-resolution simulations with gas cooling, star formation, 
and stellar feedback processes [50-53]. Importantly, ref. [53] does an extensive error analysis 
of their numerical results. In contrast, some studies' simulations have claimed that baryonic 
effects may have the opposite effect, reducing the dark matter density in the central region 
via dark matter expansion from stellar and gas feedback outflows, and producing flat or 
nearly-flat density cores at up to 2-3 kpc in size for an approximately Milky Way size halo, 
fit by cored isothermal or Burkert profiles [54-56] . It has been shown in ref. [57] that the S2 
or Plummer force softening must be a factor of ~5 times smaller than the scale of interest for 
the inner profile of dark matter halos in order to achieve greater than 5% accuracy in radial 
accelerations of particles, with an ideal time-step algorithm choice. Poor force resolution has 
been shown to generally lead to artificially lower central densities [58]. The gravitational 
softening length used in the simulations in ref. [55] is 0.5 kpc, so it is questionable to draw 
conclusions at the claimed ~kpc core scale. In the recent work of ref. [56], the gravitational 
softening length is 0.3125 kpc [59], and therefore conclusions of the inner 1-2 kpc are also 
difficult to make with confidence. Furthermore, ref. [56] finds evidence for a core to be pro- 
duced only in the more extreme feedback High Feedback Run, while the Low Feedback Run 
found adiabatic contraction that steepened the dark matter profile. 

In summary, work indicating the presence of dark matter density cores in numerical 
simulations are at the edge of their resolution limits. It is important to consider that if it 
becomes firmly established from numerical simulations that there is necessarily a constant- 
density core in the Milky Way at the > 450 pc scale, then the HESS GC limits presented 
here are not applicable due to the observation's background subtraction method. This was 
also noted by the HESS collaboration work [3]. However, at this time robust numerical 
simulations predominantly indicate contraction and steepening of the central density profile 
due to baryonic effects, and the adoption of the non- contracted NEW or Einasto profile here 
is conservative relative to the steeper profiles. 

Below we show that, in the case of a non- adiabatically-contr acted NEW or Einasto dark 
matter halo profile, the HESS GC observation provides a strong limit on the cross-section 
of high-mass WIMPs' annihilation into several Standard Model channels. Furthermore, we 
show how the dark matter annihilation interpretation of the PAMELA excess and Fermi-LAT 
e'^e~ feature signals is excluded at above 95% CL in many cases. Dark matter interpretations 
of these signals are in tension with the HESS GC observations for two-body standard model 
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particle final states and gauge-boson mediated four-lepton final states, when also including 
the constraints for Fermi-LAT observations toward local dwarf galaxies. Importantly, the 
HESS limits presented here for XDM are more constraining than the thermal relic density, 
CMB, self-interaction bounds, and naturalness bounds. 



2 Data Analysis 

The data used in this paper comes from the HESS collaboration analysis of the GC from 
ref. [3]. The events analyzed are those from 112 hours of live time from the HESS very-high 
energy gamma-ray instrument with zenith angles smaller than 30° which were within the cen- 
tral 4° of the HESS field-of-view. Contamination of the dark matter signal due to the Galactic 
plane is excluded by masking the regions with Galactic latitude |6| < 0.3° and an additional 
mask within 0.3° of the extended source HESS J1745-303 ((6,£) = (-0.6°, 358.71°)) [60]. 

The source region is defined by 0.02° x0.02° pixels that lie within 1° of the GC, do not lie 
within the mask, and have a well-defined background. The background region is determined 
by choosing a telescope pointing position within 1.5° of the GC and rotating each pixel in 
the source region (with masked pixels removed) by 90°, 180°, and 270° about the telescope 
pointing position. Any of these pixels which is further than 1° from the GC and does not 
lie within the mask is considered background. Pixels within the inner 1° for which there are 
no corresponding background pixels are excluded from the analysis. This process was then 
repeated for multiple telescope pointing positions, calculating the differential flux in both the 
background and source regions. For additional details on the HESS GC analysis, see ref. [3] . 



3 Dark Matter Annihilation Limits from the HESS GC Observations 
3.1 Gamma-Ray Emission from Annihilating Dark Matter 

A robust calculation of the expected final state radiation from dark matter annihilation 
requires accurate quantification of the dark matter source as well as the products in the 
final state gamma-ray radiation chain. The differential flux per solid angle for a dark matter 
candidate with cross-section {aAv) over a solid angle AO, is 



d£ _ (aAv) Jaq 1 dN^ 
dE~ 2 Jo 47rM2 dE 



(3.1) 



where dN^/dE is the gamma-ray spectrum per annihilation and My- is the dark matter 
particle mass. The integrated mass density squared along line-of-sight x, averaged over the 
solid angle of the observation region Ail is defined as 



where distance from the GC is given by 



d^ j dx p^{r^^i{b,l,x)) , (3.2) 



rgai(6, x) = ^Jr% - 2xRq cos{i) cos(6) + . (3.3) 

A normalization constant Jq = 1/ [8.5 kpc(0.3 GeV cm^'^)^] is chosen to make J dimension- 
less, but the final flux calculation is independent of the choice of Jq. The coordinates b and i 
are the Galactic latitude and longitude, respectively. Following the HESS Collaboration, the 
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dark matter profiles p(r) we choose are the non-adiabatically-contracted Einasto and NFW 
models of ref. [61] with the local dark matter density pq = 0.389 GeV cm~^ (0.385 GeV cm~^) 
for an NFW (Einasto) profile [62]. Since there were several telescope pointing positions not 
given explicitly in ref. [3], we cannot and did not independently calculate J ah but adopt 
those of the HESS Collaboration. Importantly, as discussed in the introduction and in 
the HESS Collaboration study, in the case of a cored-isothermal or Burkert profile with a 
constant-density core that extends at or beyond ~450 pc, then the background subtraction 
region would have an identical annihilation signal as the signal region, and no constraint can 
be placed on dark matter annihilation by this method. We have checked that the J-values 
for the single pointing shown in figure 2 of ref. [3] are approximately their J's. Averaging 
over telescope pointing positions, the HESS analysis calculated the dark matter J-values for 
the source (Jg) and background {J h) regions using the NFW and Einasto profiles [3]: 



JNFW ^ igo4 


(3.4) 
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(3.5) 


J Einasto ^ 3^43 


(3.6) 


J Einasto ^ ^535 ^ 


(3.7) 



Note that toward the GC we adopt zero astrophysical substructure boost. However, there 
may be a substructure contribution to the flux of order 4% to 0.04% in the inner ~1° view 
of the GC [63]. 

3.2 Calculation of Dark Matter Spectra 

To calculate the photon spectrum for a particular WIMP annihilation channel, we use PYTHIA 
6.4 to simulate the photon radiation of charged particles as well as decays of particles such 
as the vr'^ [64]. Specifically, we run pythia to simulate an e'^e~ collision at a center of mass 
energy of 2M^ through a Z' to a final state that corresponds to the annihilation products 
of the dark matter. For WIMP annihilation to final states through light gauge bosons (such 
as the 4e and 4;U channels) we have the Z' decay into two scalar (Hq) states, each of which 
annihilates into two leptons. For the XDM case, we employ a light gauge boson which has 
branching fraction to e~^e~, p,~^fi~, and 7r"'"7r~ set by the particular XDM model. For the 
models with light gauge bosons, photons can only be radiated at fairly low energy in the gauge 
boson rest frame. Therefore, in the center-of-mass frame of the dark matter annihilation, the 
average number of hard photons is significantly reduced in comparison to direct annihilation 
to a two-body standard model final state. 

We turn off initial state radiation such that all photons only come from the radiation 
or decay of the dark matter annihilation products. We turn on the decays of particles which 
are not decayed with the default PYTHIA settings, namely muons, charged pions, and charged 
kaons. Additionally, we turn on the muon decay channel fi~ — )• e~f^Pe7i with the standard 
branching fraction of 0.014 [65]. Using a large sample of events for each final state and each 
value of My., the number of photons in the final state in a given logarithmic energy bin is 
counted and averaged over the number of events, yielding the average number of photons in 
that energy bin per annihilation event. 

For WIMP annihilation to final states through extremely light gauge bosons (m^ < 
0.5 GeV), the pythia calculation becomes impractical due to low-energy cutoffs for radiation 
processes. For such annihilation channels (with = 0.25 GeV and m,^ = 0.35 GeV) we 
instead use the analytic formulae given in Appendix A of ref. [66]. For the 4e state, all 
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photons come from radiation off the final electrons and positrons. For the 4^ state, photons 
come from radiation off the final electrons and positrons as well as the radiation from the 
muons before they decay. We have verified the analytic spectra model by comparison to 
PYTHIA spectra for heavier gauge bosons. 

3.3 Limits on the Dark Matter Annihilation Cross-Section 




Figure 1. Shown are the data from the HESS GC observation, for which the signal and background 
are consistent with no difference in flux at the level of x^/DOF = 0.75. For comparison, two possible 
dark matter signals are shown as well. Both dark matter signals are for a 1.2 TeV WIMP annihilating 
via 0.25 GeV gauge bosons into two e+e^ pairs, for an NFW halo profile. The red (lower) histogram 
is the signal expected for a dark matter cross-section {(Jav) = 9 x 10~^^ cm'^ s~^, which has a total 
Ax^ = 2.79. The magenta (upper) histogram is the signal expected for a dark matter cross-section 
((Taw) = 5 X 10^^^ cm'^ s~^ and mass — 1.2 TeV, with a Ax^ — 41.2. These signals correspond 
with the upper end of the red Fermi-LAT plus PAMELA region, and lower end of the purple rectangle 
in figure 4(a), in the case of assuming a subdominant substructure contribution, and are excluded at 
greater than 95% CL. Note that the HESS GC observation extends to higher energies than shown 
here. 

Figure 3 of ref. [3] shows the observed source and background fluxes in 35 energy bins 
from 0.28 to 31 TeV. The two regions are consistent with each other, with zero difference in 
flux having a x^/DOF = 0.75. Therefore, any dark matter signal must be small enough that 
the source region does not have appreciably greater gamma-ray flux, within errors, than the 
background region. Figure 1 illustrates the difference between the lack of a dark matter signal 
from HESS and two representative dark matter signals. Using the dark matter J-values from 
eqs. 3.4-3.7, we have derived 95% one-sided confidence-level (CL) (corresponding with a 90% 
two-sided CL) constraints on the dark matter cross-section as a function of mass for several 
key annihilation channels, using the total iii bins. 

In the figures, the light blue cross-hatched region is excluded at 95% CL for both the 
NFW and Einasto dark matter profiles. The singly-hatched light blue regions are excluded 
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at 95% CL for the Einasto profile but not for the more conservative NFW profile. The purple 
regions are the 95% CL limits from a combined analysis of ten dwarf spheroidal galaxies from 
the Fermi-LAT Collaboration [15]. In particular, note that in figure 2(a) and figure 3(b) we 
have extended the mass of the dark matter down to 10 GeV in order to show the exclusion 
of a standard thermal relic below 27 GeV (37 GeV) for the bb (t'^t~) annihilation channel 
by the Fermi-LAT stacked analysis of dwarf galaxies [15]. The dark pink regions represent 
the annihilation cross-section for a thermal relic, {(tav) ~ 3 x 10~^^ cm'^ s~^. The dark 
matter constraints for the standard WIMP annihilation channels bb, tt, and W~^W~ are 
shown in figure 2. We show the tt channel since this should dominate for high mass dark 
matter (M^ > 200 GeV). Our bb constraints are consistent with the quark channel limits 
from HESS [3]. As a comparison to the thermal cross-section, figure 2(c) also includes the 
expected cross-section for a non-thermal wino-like neutralino (the dashed red line) [67] . 

Figure 3 shows the dark matter limits for dark matter annihilating directly into leptons, 
whereas in figure 4, the dark matter annihilates into two gauge bosons cp of mass m,^ = 
0.25 GeV which then decay into leptons. (Note that the 4e case requires an ad hoc requirement 
of gauge boson decay into electrons and not muons. Neither the constraints nor signals 
would be significantly different in the case of, e.g., = 0.2 GeV, where annihilation to 
e'^e~ would be energetically required.) Heavy dark matter masses annihilating primarily 
into leptons are particularly interesting in the context of the PAMELA positron excess [16] 
and the e~^e~ feature seen by the Fermi-LAT [20]. Such leptonic annihilation channels with 
a cross-section much larger than that of a thermal relic have been studied as the source of 
these anomalies [10, 22]. In figures 3 and 4 we have included one dark matter annihilation 
interpretation of the PAMELA excess in light pink and the analogous interpretation of the 
Fermi-LAT feature in red [10]. Recent analysis of the positron fraction by the Fermi-LAT 
Collaboration consistent with the PAMELA excess continues to rise up to energies of 180 
GeV bin center energy, ruling out PAMELA regions below ~ 160 GeV from being 
consistent with dark matter annihilation [23]. In figures 3(a) and 4(b) we include fits to the 
interpretation of the PAMELA excess in light green outline and the Fermi-LAT e~^e~ feature 
in dark green outline from ref. [22]. In panel 4(b), for the NFW profile case, we exemplify 
the strength of the limit by plotting the 95%, 99.7%, and 99.9999% CL hmits as dashed, 
dot-dashed and solid lines, respectively. 

In figure 5 three benchmark XDM models which are consistent with a combination of 
the PAMELA signal and the Fermi-LAT feature are shown, with red representing the 68% 
CL region and light pink representing the 95% CL region. Figure 5(a) contains the regions for 
annihilations which go 50% into e~^e~ and 50% into ^jl'^ through two intermediate gauge 
bosons of mass = 0.35 GeV; figure 5(b) contains the regions for annihilations which go 
33% into e'^e~ , 33% into IJ''^H~ , and 33% into vr+vr" through two intermediate gauge bosons 
of mass 7710 = 0.58 GeV; and figure 5(c) contains the regions for annihilations which go 25% 
into e+e", 25% into fJ''^fi~, and 50% into tt+tt" through two intermediate gauge bosons of 
mass 7770 = 0.90 GeV. 

When comparing to signal-fit regions from other work, a scaling of Pq/Pq has been 
done to normalize that work's local dark matter density po to the one adopted here. In 
the literature, signals for the PAMELA and Fermi-LAT excesses are ascribed boost factors 
for dark matter annihilation that can include both astrophysical substructure boosts and 
particle physics boosts, as well as the enhancement of latter due to the former. The local 
boost can be separated as i?iocai = BgBp into the substructure boost Bs and the particle 
boost Bp (when ignoring the enhancement of the latter to to the former). To be clear and 
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conservative, we employ a relatively strong local substructure boost, Bs = 1.57, using that 
expected from unresolved substructure calibrated to the Via Lactea II simulations [38, 69]. 
We incorporate the fit regions from ref. [22] such that 

^''=( o.3GeVcm-0 '^-^-- P'**' 
Both -Biocai ^-nd Ep designate the scaling factor from (aA'v) = 3 x 10~^® cm^ s~^, as shown 
in our figures. 

Relative boosts between the GC and local effects can alter the relative constraints 
between the annihilation in the GC versus local cosmic ray signals in Sommerfeld-enhanced 
annihilation models, where Sommerfeld-enhancement is greater in low- velocity substructures. 
As discussed in ref. [37], the boost toward the GC relative to the local is 

BqC _ Sy^r=0)/Svr^l50 km/s 

-Biocal 1 + {Sv^o/Sy^i5o km/s)A(8.5 kpc) 

where S is the Sommerfeld-enhancement in a given model at a distance r from the GC, for 
a given velocity dispersion of the dark matter v. Here, A = 1 — Bg, which we take to be 
A(8.5 kpc = 0.57. 

For specificity, we adopt m,^ = 0.25 GeV, which corresponds to Sv^o/ S^r^i^Q km/s ~ 5, 
from figure 1 of ref. [37]. This model is a 4e channel case shown in figure 4(a). The purple 
rectangle shows the range of annihilation cross sections for a 1.2 TeV Sommerfeld-enhanced 
scenario consistent with the thermal relic density, the CMB, self-interaction bounds, and 
naturalness, for A from 10~^ to unity, corresponding to -Biocal from 30 to 300 [37]. The 
scaling of the signal toward the GC relative to the local boost is set explicitly by the local 
substructure boost A, and the value of S'^^^.^o)) which is determined by the velocity dispersion 
of dark matter in the ~450 pc of the HESS observation of the GC (and not exactly at r = 0). 
It has been shown that the dark matter velocity rapidly decreases toward the GC in Milky 
Way scale halos [70], and the limit t; — )• is potentially appropriate for 5^,(^=0), which 
we adopt in one case of our limits. In this GC low velocity limit case, one can take the 
substructure dominant versus subdominant cases: 

(A = 10-4), , , 

^ ^' 3.10) 

(A = l). 

Baryonic effects have been found to enhance the velocity of the dark matter toward the 
GC to make it comparable or greater than to that at the solar distance [50, 71], such that 
Sy(^r=o) ~ 'S',;~i5o km/s [72], and in that scenario the substructure limiting cases are: 

Bgc{v ~ 150 km s-i) ^ [ 1 (A = 10-4), 




-Blocal [0.17 (A = l) 



We designate these ranges of constraints as bars and arrows for the corresponding four cases 
in figure 4(a), with the f; — )• GC in green and v ~ 150 km GC in dark blue. 

Ref. [37] claims that A(8.5 kpc)/A(r = 0) ~ 20 is relevant for the comparative con- 
straints between the GC and the local boost, which strictly is only the case for strong local 
substructure domination and high velocities in the GC. More generally, it is eq. (3.9) and 
the local substructure boost that sets the scaling, and is what we adopt. Note that any 
substructure in the dark matter of the GC region would enhance the constraints here. We 
ignore this enhancement to be conservative, i.e., we only include the annihilation from the 
smooth component in the GC region. 



-8- 



3.4 Discussion 

For the bb annihilation channel (figure 2(a)), the HESS GC constraints are stronger than 
the constraints from the Fermi-LAT analysis of dwarf spheroidal galaxies for > 900 GeV 
when adopting a non-cored Einasto or NFW dark matter profile. The limits on this channel 
for an Einasto dark matter profile are within an order of magnitude of the thermal cross- 
section for 2 TeV < M^. < 5 TeV. Similarly, the W~^W~ annihilation channel (figure 2(c)) 
has stronger constraints than the Fermi-LAT dwarfs for > 800 GeV and limits {(Jav) < 
3 X 10~^^ cm^ s~^ for 1 TeV < < 6 TeV. The ti annihilation channel (figure 2(b)) has 
somewhat weaker constraints, limiting (cta'u) < 5 x 10^^^ cm^ s~^ for 2 TeV < < 10 TeV. 

The light pink PAMELA excess region of ref. [10] and the light green PAMELA excess 
region of ref. [22] are both excluded above ~ 400 GeV by the HESS GC data for the 
fJ-^f^~ annihilation channel (figure 3(a)) when adopting a non-cored Einasto or NFW halo 
model. Also for the ^^/U^ channel, the red Fermi-LAT feature region of ref. [10] and the 
dark green Fermi-LAT feature region of ref. [22] are excluded by the HESS GC data when 
adopting either halo model. 

In the r+T~ annihilation channel (figure 2(b)), the HESS GC observation excludes a 
cross-section of {(tav) > 4 x 10^^^ cm^ s~^ for a dark matter mass ~ 1 TeV, when 
adopting an Einasto profile, within a factor of ~2 of the thermal cross-section. The HESS 
GC excludes the light pink PAMELA excess region above My. ~ 400 GeV and the Fermi- 
LAT dwarf analysis excludes the r+r" channel below ~ 400 GeV, so this model for 
the PAMELA excess has been ruled out at all dark matter masses when adopting an NFW 
profile. The red Fermi-LAT feature region is also excluded, at greater than 99.9999% CL, 
in the t~^t~ channel for both the Einasto and NFW dark matter profiles, consistent with 
previous results [7, 10, 12]. 

For the 4e annihilation channel (figure 4(a)), the HESS GC observation constrains the 
light pink PAMELA excess region of ref. [10] above ~ 600 GeV for the Einasto halo 
model, and is constrained above ~ 500 GeV for the NFW halo model. Similarly, the 4/x 
annihilation channel (figure 4(b)) has both the light pink region and the light green PAMELA 
excess region of ref. [22] constrained above My. ~ 1.5 TeV for the Einasto halo model, and 
above w 3 TeV for the NFW halo model. 

In the 4e channel (figure 4(a)), the purple rectangle shows the range of boost factors for a 
1.2 TeV Sommerfeld-enhanced model consistent with the thermal relic density, the CMB, self- 
interaction bounds, and naturalness, for A conservatively from 10^^ to unity [37]. Note that 
the gamma-ray annihilation constraints from the HESS GC on the annihilation cross-section 
and boost factor, when adopting an NFW or Einasto profile, is often stronger than those from 
the thermal relic density, the CMB, self-interaction bounds, and naturalness considerations. 
Specifically, the limits to boosts toward the GC are Bqq ^ 27, which for the comparable 
central velocity case excludes models with local substructure A < 0.025 (0.17) for .Biocai = 
30 (50). Our canonical case for A w 0.57, with GC dark matter velocities comparable to 
that locally, is unconstrained at 95% CL. Note that any substructure within the inner degree 
of the HESS GC observation would further enhance these limits due to the corresponding 
Sommerfeld-enhancement saturation in the GC substructure. The substructure boost in the 
inner degree has been seen in simulations at the level of A(r ~ 0.15 kpc) ~ 4 x 10~^ to 4 x 
10"^ [63]. 

The XDM annihilation models: Br(e'''e~) = Br(/i'^^~) = 0.5, niif, = 0.35 GeV (fig- 
ure 5(a)); Br(e+e-) = Br(^+/i~) = Br(7r+7r-) = 0.33, = 0.58 GeV (figure 5(b)) ; and, 
Br(e+e-) = Br(;U+^-) = 0.25,Br(7r+7r-) = 0.5, = 0.90 GeV (figure 5(c)) are excluded 
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at greater than 95% CL when adopting an Einasto halo profile. However, when adopting 
an NFW halo profile, these XDM models have a region from 0.8 TeV < < 1 TeV below 
{(JAv) ~ 1.8 X 10~^^ cm^ s~^ which remains consistent at 95% CL. 

4 Conclusions 

The HESS telescope's observations toward the Galactic center present the strongest con- 
straints on WIMP dark matter annihilation into Standard Model particles for > 900 GeV, 
given a non-adiabatically-contracted NFW or Einasto profile for the Milky Way dark matter 
profile. As discussed in the introduction, the HESS GC observation is not sensitive to dark 
matter annihilation for large constant-density cored dark matter profiles. If such profiles are 
established to be valid in the Milky Way, the HESS GC observation provides no empirical 
constraint in these cases given the background subtraction method. If adiabatic contraction 
and steepening of the dark matter density profile is well established, then the constraints 
could become stronger. We find that constraints on annihilation into bb final states comes 
within an order of magnitude of the canonical thermal cross section, while tt is within a factor 
of ~20. The T^T~ channel is within a factor of ~2 of the thermal cross-section. This bodes 
well for the future Cerenkov Telescope Array's potential impact at constraining canonical 
thermal WIMP dark matter this mass scale [73]. 

We also examine constraints on Sommerfeld-enhanced dark matter annihilation models 
which produce the PAMELA positron excess and Fermi-LAT e~^e~ spectral feature, includ- 
ing XDM. The models with pure leptonic modes t~^t~ and 4e are excluded at greater than 
95% CL by the HESS GC when adopting NFW or Einasto dark matter halo profiles. For 
other cases, (e.g. 4/i, 2e2/i , e~ fi'^ fi^ tt^ ir^ ) , the models are in tension with the HESS 
GC observations, with portions of the signal 95% CL parameter space are excluded at the 
95% CL level. Significantly, the exclusions presented here from HESS GC gamma-ray obser- 
vations on Sommerfeld-enhanced upper boosts are more constraining, in many cases, than 
prior constraints from diffuse gamma-rays [38, 39], relic density considerations [28-32], the 
CMB [29, 33-35], halo shapes [30, 31], and naturalness [36, 37]. In the case where the veloc- 
ity dispersion of dark matter in the center of the Galaxy is comparable to that locally, the 
supplemental Sommerfeld-enhancement from local substructure and not at the GC weakens 
the constraints in certain models. 

Gamma-ray astronomy has produced the most stringent constraints on the canonical 
thermal WIMP model's annihilation cross-section, with Fermi-LAT's stacked observations of 
dwarf galaxies being the most constraining at low masses, and HESS's observations of the 
Galactic center being the most constraining at higher masses. With further observation and 
new technologies, the nature of dark matter may be revealed by gamma-ray telescopes. 
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Figure 2. Shown are the constraints on dark matter in three canonical annihilation channels: (a) bb; 

(b) tt; (c) W^W^. The regions are labeled according to their constraining observations as described 
in the text: "HESS GC" are the 95% CL limits from the HESS analysis of the GC. The double hatched 
region is constrained for both the Einasto and NFW halo models, and the single hatched region is 
constrained for only the Einasto halo model. The regions labeled "Fermi Dwarfs" are the 95% CL 
limits from the Fermi-LAT collaboration analysis of dwarf spheroidals. In the W~^W~ channel, panel 

(c) , the mass for a non-thermal wino-like neutralino is shown as a thick-dashed red line [67]. For 
comparison, we plot the Scr limits from ref. [68] for their analysis of the Fermi-LAT observation of 
the 3° X 3° region around the Galactic Center as dashed (black) lines in all panels for the respective 
channels. 
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Figure 3. Shown are the constraints on dark matter in two leptonic annihilation channels: (a) fi'^ fi ; 
(b) r+T^. The regions are labeled according to their constraining observations as described in the 
text: "HESS GC" are the 95% CL limits from the HESS analysis of the GC. The double hatched 
region is constrained for both the Einasto and NEW halo models, and the single hatched region is 
constrained for only the Einasto halo model. The regions labeled "Fermi Dwarfs" are the 95% GL 
limits from the Fermi-LAT collaboration analysis of dwarf spheroidals. The light pink shaded region 
is consistent with a dark matter interpretation of the PAMELA signal and the dark red shaded region 
is that from the Fermi-LAT e+e~ feature from ref. [10]. In the channel, panel (a), the light 

green outlined region is consistent with a dark matter interpretation of the PAMELA signal and the 
dark green outlined region with that of the Fermi-LAT e~^e~ feature from ref. [22]. PAMELA regions 
below = 160 GeV are ruled out by the rise in the positron fraction seen by the Fermi-LAT [23]. 
In panel (b), to illustrate the strength of the HESS GG limits, we show for the NFW profile the 95%, 
99.7% and 99.9999% GL limits in dotted, dot-dashed and solid lines, respectively. For comparison, 
we plot the 3a limits from ref. [68] for their analysis of prompt and inverse-Gompton emission in the 
Fermi-LAT observation of the 3° x 3° region around the Galactic Genter as dashed (black) lines in 
both panels for the respective channels. 
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Figure 4. Shown are the constraints on dark matter in two annihilation channels: (a) annihilation 
into two e+e^ pairs via two intermediate 0.25 GeV gauge bosons </) (note that this case requires 
an ad hoc requirement of decay into electrons and not muons); (b) annihilation into two /x"*" and 
two /i~ via two intermediate 0.25 GeV gauge bosons 4>. The regions are labeled according to their 
constraining observations as described in the text: "HESS GC" are the 95% CL limits from the HESS 
analysis of the GC. The double hatched region is constrained for both the Einasto and NEW halo 
models, and the single hatched region is constrained for only the Einasto halo model. The light pink 
shaded region is consistent with a dark matter interpretation of the PAMELA signal and the dark 
red shaded region is that of the Fermi-LAT e~^e^ feature from ref. [10]. In the 4e channel, panel (a), 
the purple rectangle demonstrates the range of Sommerfeld-enhanced cross-sections consistent with 
constraints from thermal relic density, the CMB, self-interaction bounds, and naturalness [37]. The 
green (blue) bars and arrows show the HESS GC limits for the two cases of velocity dispersions of 
w — >■ (w ^ 150 km s~^), respectively, with the upper and lower bars for each color representing the 
two local substructure boost limits, as described in the text. In the 4/j, channel, panel (b), the light 
green outlined region is consistent with a dark matter interpretation of the PAMELA signal and the 
dark green outlined region is that for Fermi-LAT e+e" feature from ref. [22]. PAMELA regions below 
« 160 GeV are ruled out by the rise in the positron fraction seen by the Fermi-LAT [23]. 
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Figure 5. Shown are the constraints on dark matter in three XDM annihilation cases: (a) annihilation 
into 50% e+e~ and 50% via two intermediate 0.35 GeV gauge bosons 0; (b) annihilation into 

33% e+e~, 33% fJ.~^fJ.~ , and 33% tt+tt" via two intermediate 0.58 GeV gauge bosons 0; (c) annihilation 
into 25% e+e~, 25% and 50% tt+tt^ via two intermediate 0.90 GeV gauge bosons cj). The 

regions are labeled according to their constraining observations as described in the text: "HESS GC" 
are the 95% CL limits from the HESS analysis of the GC. The double hatched region is constrained 
for both the Einasto and NEW halo models, and the single hatched region is constrained for only the 
Einasto halo model. The light pink and dark red shaded regions are the 68% and 95% CL regions 
consistent with the XDM dark matter interpretation of a combination of the PAMELA signal and the 
Eermi-LAT e+e^ feature from ref. [36]. 
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